Moments of pseudoscalar meson distribution amplitudes from the lattice 
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Based on lattice simulations with two flavours of dynamical, 0(a)-improved Wilson 
fcrmions we present results for the first two moments of the distribution amplitudes of 
pseudoscalar mesons at several values of the valence quark masses. By extrapolating our 
results to the physical masses of up/down and strange quarks, we find the first two moments 
of the K + distribution amplitude and the second moment of the ir + distribution amplitude. 
We use nonperturbatively determined renormalisation coefficients to obtain results in the 
MS scheme. At a scale of 4 GeV 2 we find = 0.201(114) for the second Gegenbauer 
' moment of the pion's distribution amplitude, while for the kaon, af = 0.0453(9) (29) and 

af = 0.175(18)(47). 
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I. INTRODUCTION 



Qh. In recent years exclusive reactions with identified hadrons in the final and/or initial state are 

attracting increasing attention 0] . The reason for this interest is due to the fact that they are dom- 
inated by rare configurations of the hadrons' constituents: either only valence-quark configurations 
contribute and all quarks have small transverse separation (hard mechanism) [2J, y, |j, la, la, LD, |8|, |{| , 
or one of the partons carries most of the hadron momentum (soft or Feynman mechanism). In 
both cases, the information about hadron structure is new and complementary to that in usual 
inclusive reactions, the prominent example being the deep-inelastic lepton hadron scattering. 

Hard contributions are simpler to treat than their soft counterparts and their structure is well 
understood, see e.g. Ref. ^fl*] for a recent discussion. They can be calculated in terms of the hadron 
distribution amplitudes (DAs) which describe the momentum-fraction distribution of partons at 
zero transverse separation in a particular Fock state, with a fixed number of constituents. DAs are 
ordered by increasing twist; the leading twist- 2 meson DA, <pu, which describes the momentum 
distribution of the valence quarks in the meson II, is related to the meson's Bethe-Salpeter wave 
function <pn,BS by an integral over transverse momenta: 

Km 



/|K_L|<jU 
d 2 k± 4>u,bs(x, k±) 



Here x is the quark momentum fraction, Z<i is the renormalisation factor (in the light-cone gauge) 
for the quark-field operators in the wave function, and \x denotes the renormalisation scale. In 
particular the leading-twist DA of the pion and of the nucleon have attracted much attention in 
the literature. Furthermore, SU(3) flavour symmetry breaking effects in the DAs of strange mesons 
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are important for predictions of the exclusive B-decay rates (e.g. B — > K, K*) in the framework of 
QCD factorisation HH, perturbative QCD 0, soft-collinear effective theory (SCET) HE El or 
light-cone sum rules, e.g. 3, IE) 01- I n some cases, for instance weak radiative decays, B — > pj 
vs. B — > K*j, the uncertainty in SU(3) breaking is actually the dominant source of theoretical 
error. 

The theoretical description of DAs is based on their representation |E IE IE IE IE 0, IE IE] as 
matrix elements of a suitable nonlocal light-cone operator. For example, for positively charged 
pions or kaons one defines 

(0\q(-z)^ 75 [-z,z}u(z)\U + (p)) =if uP ^ j'^e-^M^^ 2 ), U) 

where q = d,s, z^ is a light-like vector, z 2 = 0, \—z, z] is the straight-line-ordered Wilson line 
connecting the quark and the antiquark fields and /n is the usual decay constant f w = 132 MeV, 
fx = 160 MeV [18(| . The physical interpretation of the variable £ is that x = (1 + £)/2 and 
1 — x = (1 — £)/2 are the fractions of the meson momentum carried by the quark and antiquark, 
respectively. The definition in Q implies the normalisation 

J d£0n(£,M 2 ) = l. (2) 

For brevity, below we often drop the subscript and write <fi instead of 4>u unless we are referring to 
a specific meson. 

A convenient tool to study DAs is provided by the conformal expansion 0, EE El- The 



underlying idea is similar to the partial-wave decomposition in quantum mechanics and allows one 
to separate transverse and longitudinal variables in the Bethe-Salpeter wave-function. The depen- 
dence on transverse coordinates is formulated as a scale dependence of the relevant operators and 
is governed by renormalisation- group equations. The dependence on the longitudinal momentum 
fractions is described in terms of Gegenbauer polynomials C n (£) which are nothing but irreducible 
representations of the corresponding symmet ry g roup, the collinear conformal group SL(2,M). 



aresentations ol tne corresponding symmetry 
In this way one obtains [j, [E 0, 0, 3> iE 2. 



71=1 



(3) 



To leading- logarithmic accuracy (LO), the (non-perturbative) Gegenbauer moments a n renormalise 
multiplicatively with 

a n (^) = L^^Ma n (ri), (4) 
where L = a s (p 2 )/a s (fj,Q), 0o = 11 — 2iVy/3, and the anomalous dimensions 7^ are given by 

^ 8 c F (gi-^- 2(n+i ; (n+2) ) (s, 

with Cp = 4/3. Note that the multiplicative renormalisability in leading order is not an acci- 
dent: It relies on the fact that the tree-level counterterms retain the symmetry properties of the 
Lagrangian |2l| . 

Since the anomalous dimensions increase with spin, the higher-order contributions in the Gegen- 
bauer expansion are suppressed at large scales so that asymptotically only the leading term survives 

<M£,/i 2 ^oo) = <M£) = !(l-£ 2 ). (6) 
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For this reason, one usually assumes that the conformal expansion is well convergent at all scales 
of practical interest, and retaining the first few terms only in the conformal expansion provides one 
with a reasonable approximation for convolution integrals of the type d£ 0(£)/(l — £) that one 
encounters in many applications. 

To next-to-leading order (NLO) accuracy, the scale dependence of the Gegenbauer moments is 
more complicated and reads 0, [ij, EF 



/ 2\ 71—2 
NLO , a s\^ ) S^„_f..2\ pNLO 



Air 



(7) 



k=0 



Here we adopt the usual convention that an empty sum is equal to zero. Moreover, ao = 1 and 



E 



NLO 



i + 



7^0 - 7^i 



8vr/3 n 2 



a s (n 2 ) -ot s {nl) 



(8) 



where 7^ are the diagonal two- loop anomalous dimensions [26j], f3\ = 

mixing coefficients d^, k < n — 2, are given in closed form in Ref. 24], 
Ref. 27] for a recent compilation. For the lowest moments 



102 - (38/3) N f , and the 
2II . see also, for instance, 



0, 



and 



''20 



7i 
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~ 243 ~ "8T 



N f 



(9) 



(10) 



If the mass difference between the u and d quarks is neglected, G-parity implies that the pion DA 
^7r(£)M) is an even function of £, i.e all odd moments in ^, a^ra+i' van ish. The K-meson DA need 
not be even, and the calculation of af will be one of our goals. 

The coemcients a n at some reference scale //o are nonperturbative quantities and have to be 
evaluated using a nonperturbative technique or must be extracted from experiment. In historic 
perspective, most of the discussion over the years was centered on a particular model of the pion DA 
proposed by Chernyak and Zhitnitsky in 1982 on the basis of their calculation of using QCD 
sum rules [28j. Using this model and assuming dominance of the hard rescattering mechanism 
in exclusive reactions, Chernyak and Zhitnitsky were able to describe an impressive amount of 
experimental data that were available at that time [2^. 

Since then, the original argumentation by Chernyak and Zhitnitsky and the model itself have 
been largely discredited. Three different approaches have been used: direct calculations using 
QCD sum rules, pioneered in 12811 : analysis of experimental data on the pion electromagnetic and 
transition form factors (e.g. 30]) and the B weak decay form factor, using light-cone sum rules; 
and lattice calculations. The summary of these results is presented in Table 2 of Ref. 0|; see 
also, for instance, Refs. [27], E3] for another recent compilation. It turns out that can only be 
determined with large errors, whatever approach is chosen. A fair average is probably 



4(4 GeV 2 ) =0.17 ±0.15. 



(11) 



The if -meson DA has attracted comparatively less attention. The numerical value of the first 
moment af was the subject of significant controvers y u ntil rece ntly . The existing estimates are all 
obtained using different versions of QCD sum rules 3^, 3 3^, 3(| and yield an average ET 



,A 



(4 GeV 2 



0.05 ±0.03. 



(12) 
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For the second moment, the old estimate by Chernyak and Zhitnitsky was ja\ = 0.59 ± 0.04. 
Two recent sum rnle calculations 0,0, including radiative corrections to the sum rules, give 
however af/aj — 1 pointing towards a very small SU(3) violation in this coefficient. 

Estimates of yet higher-order Gegenbauer coefficients are rather uncertain. The light-cone sum- 
rule calculations of the transition form factor F W11 * in Refs. 0, 37, H, 3l| suggest a negative value 



for aj, which is consistent with the result aJ(lGeV 2 ) > —0.07 obtained in Ref. |40J. However, 
this conclusion may be premature because yet higher moments have been omitted (however, in 
Ref. 41] they are estimated to be small). Moreover, there does not seem to be any convincing 
method to estimate the uncertainty due to the model dependence of the analysis. While it seems 
that distribution amplitude moments beyond the second are extremely difficult to access on the 
lattice, it might be possible to estimate them using the transverse lattice approach [42| or from the 
operator product expansion in lattice QCD 0,0]. 



Last but not least, we have to mention the estimate of the pion DA in the middle point where 
the momentum is shared equally between the quark and the antiquark 45] 

^(f = 0) = 0.6 ±0.15 , (13) 

and the measurement of the pion DA in diffractive dijet production by E791 |46j. Unfortunately, 



it turns out that collinear factorisation is broken for dijet production 



DY J 




48| , so that the interpre- 



tation of this beautiful experiment is not straightforward, see also 

The lattice computation of DAs of pseudoscalar tt and K mesons will be at the focus of this 
paper. On the lattice, we cannot compute non-local matrix elements of the form However, via 
the light-cone operator product expansion (OPE), moments (£ n ) of the DAs defined by 

r>nV) = J_da n M^» 2 ) , (14) 
are related to matrix elements of the local operators 

Off ...^(0) = i n q(0)w B w ...D„ n n(0) , (15) 

by 

<0|O{J ... M (0)|n(p)) = ifuP{^ ■ . .p M (On ■ (16) 

Here M refers to the fact that the operator is defined in Minkowski space, is the covariant 

derivative, D=D — D and {. . .} denotes the symmetrisation of all indices and the subtraction of 
traces. The moments (£ n ) are related to the Gegenbauer moments a n by simple algebraic relations: 



01 = 1(0, *2 = ^ {m 2 ) - 1) , etc. (17) 

Although the first lattice calculation of {i 2 )^ appeared almost 20 years ago 0, 0], there 
has been surprisingly little activity in this area in recent times (Hil . I54I l55| to complement other 
theoretical investigations. Our preliminary result for (£ ) m was P ^se~ed7n and we found in 
the MS scheme at fi 2 = 5 GeV 2 , {(,^ s (n 2 = 5 GeV 2 ) = 0.281(28). This represents the most recent 
lattice result. The authors of Ref. [55j, on the other hand, performed a simulation in quenched QCD 
and renormalised their results perturbatively to the MS scheme at /i 2 = 7.1289 GeV 2 , (C 2 )^ S (a* 2 = 
7.1289 GeV 2 ) = 0.280(49)1°;^° . 

This paper is organised as follows. In Sec. [H] we describe our lattice technology including the op- 
erators and renormalisation prescriptions used. Section ITTT1 contains our numerical results together 
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with appropriate extrapolations towards the physical limits. Finally, in Sec. HVI we summarise our 
findings by considering the results in terms of Gegenbauer moments. Here we also discuss the 
insights that we gain on the shape of the Pion and Kaon distribution amplitudes. We tabulate our 
results in Appendix lAl 

II. LATTICE TECHNIQUES 



We define a meson two-point correlation function in Euclidean space as 
C°(t,p) = y^e-^(o {w) ... Mn} (x,i)J(0,0)t\ , 



X 



A 



^(O|0 {w ,-..Mn}(°)l n (P)> 



e -Et + ToTje -E(L t -t) 



< t < L t , 



where O^ o fln y is the Euclidean lattice transcription of Eq. (fT5|) . 



(18) 



(19) 



A = J(0)^|0), E = y "i-n + p 2 , L t is the temporal extent of the lattice and we use J(x) 

1 1 ; :/• : = q(x)^§u{x) or J(x) = A&{x) = O4 = q{x)^^^u{x) as the interpolating operator for the 
pseudoscalar mesons. The r factor tells us how the operator behaves under time reversal, t — ► Lf—t. 



We find that for tj, tjj 



while t a 4 



+. 



To increase the overlap of our interpolating operators with the ground state, we perform Jacobi 
smearing at the source |57| . while the operators we use at the sink are local. Finally, we note that 
when working with operators involving derivatives, we perform the Fourier transform in Eq. ()18|) 
at the "centre-of-mass" of the operator 



A. Choice of operators 

We need to choose the lattice operators to perform the matching of the appropriate representa- 
tions of the H(4) group — the group of Euclidean lattices relevant for our numerical computations 
— to the corresponding representations of the 0(4) group — the group of rotations and reflections 
in four Euclidean dimensions. 

For the first moment of pseudoscalar mesons containing non-degenerate mass quarks we consider 
two types of operators which we denote generically by 0^ u (fi 7^ v) and C^„, e.g. 

O a 41 = O m , (20) 
O b u = O m - 1 (o {11} + O m + {33} ) • (21) 

The first operator, 0^, requires a nonzero momentum component in the 1-direction which we 
choose as small as possible, i.e., we take p = (p, 0, 0), where p = 2tt/L s and L s is the spatial extent 
of our lattice. Using rotational symmetry, we average over the momentum choices p = (0,p, 0) and 
p = (0,0, p), using the operators in Eq. ()20|) with {41} replaced with {42} and {43}, respectively. 
The second operator, O b , can be evaluated at p = (0, 0, 0). 

In this situation, there will be no mixing with operators of equal or lower dimensions, however 
there are improvement terms that could be included _59]. Unfortunately the improvement coeffi- 
cients are not known, so we are forced to neglect their contribution, however they are expected to 
be small and as such are unlikely to affect our results. 
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For the case of the second moment, which appears for mesons with both degenerate and non- 



degenerate mass quarks, we also have two classes of operators 0®„ , O b ul/u |2fjj], e.g 



C412 = C{412} , (22) 
g£u = {m - °^ + 2 °^ . (23) 

From Eq. (|16|). we see that 0^ vp requires two non- vanishing spatial components of momentum, 
P = (PjPjO), while O b v needs only one, p = (p, 0,0) 72]. Consideration of this fact alone would 
lead one to choose O puu , since momentum components in different directions on the lattice lead to 
a poorer signal. However, lattice operators with two or more covariant derivatives can mix with 
operators of the same or lower dimension. It turns out that for forward matrix elements, O pvu 
suffers from such mixings while O pup does not. 

For matrix elements involving a momentum transfer between the two states, i.e., nonforward 
matrix elements, both operators W and O h pvv can mix with operators involving external ordinary 
derivatives, i.e. operators of the form d p d u ■ ■ ■ (q ■ ■ ■ q). For example, 0% 12 in Eq. (|2*2*|) mixes only 
with the following operator [g^J 

0«-f = d{ , dl (q 72}l5 q) . (24) 

The situation for 0\\ uv is a lot worse as it can potentially mix with up to seven different opera- 
tors 60]. While six of these operators may vanish in the continuum limit, there exists a mixing 
operator of lower dimension, and as such its contribution must be correctly taken into account 
non-perturbatively. Thus O pup offers the best possibility to extract a value of (£ 2 ) from a lattice 
simulation. 



B. Set of gauge fields 

Our gauge field configurations have been generated with two flavours of dynamical fermions, 
Nf = 2, using the Wilson gluon action and nonperturbatively 0{a) improved Wilson fermions. 
For four different values f5 = 5.20, 5.25, 5.29, 5.40 and up to four different k values per (3 we 
have generated 0(2000 — 8000) trajectories. Lattice spacings and spatial volumes vary between 
0.075-0.123 fm and (1.5-2.2 fm) 3 , respectively. A summary of the parameter space spanned by our 
dynamical configurations can be found in Table [IJ We set the scale via the force parameter, with 
ro = 0.467 fm E3|- -^ or more details regarding our definitions and conventions, see Ref. |63l ]. 

Correlation functions are calculated on configurations taken at a distance of 10 trajectories 
using 4 different locations of the fermion source. We use binning to obtain an effective distance 
of 20 trajectories. The size of the bins has little effect on the error, which indicates residual 
auto-correlations are small. 

Concerning the influence of the finite size of our lattices, our experience with other observables 
suggests that it is not significant for the ensembles considered here. However, in our simulations 
on smaller lattices (not included in the present analysis) finite size effects are to be expected, the 
study of which is under investigation. 
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TABLE I: Lattice parameters: Gauge coupling /?, sea quark hopping parameter K soa , lattice volume, the 
force scale, ro, and pion mass. The latter three are given in lattice units. 






^sca 


Volume 


r /a 


am n 


5.20 


0.13420 


16 3 x 32 


4.077(70) 


0.5847(12) 


5.20 


0.13500 


16 3 x 32 


4.754(45) 


0.4148(13) 


5.20 


0.13550 


16 3 x 32 


5.041(53) 


0.2907(15) 


5.25 


0.13460 


16 3 x 32 


4.737(50) 


0.4932(10) 


5.25 


0.13520 


16 3 x 32 


5.138(55) 


0.3821(13) 


5.25 


0.13575 


24 3 x 48 


5.532(40) 


0.25556(55) 


5.29 


0.13400 


16 3 x 32 


4.813(82) 


0.5767(11) 


5.29 


0.13500 


16 3 x 32 


5.227(75) 


0.42057(92) 


5.29 


0.13550 


24 3 x 48 


5.566(64) 


0.32696(64) 


5.29 


0.13590 


24 3 x 48 


5.840(70) 


0.23956(71) 


5.40 


0.13500 


24 3 x 48 


6.092(67) 


0.40301(43) 


5.40 


0.13560 


24 3 x 48 


6.381(53) 


0.31232(67) 


5.40 


0.13610 


24 3 x 48 


6.714(64) 


0.22081(72) 



C. Extracting the matrix elements 

We calculate the average of matrix elements computed with the following choices of momenta 

O a 



p = 


(p,0,0) 


p = 


(0,p,0) 


p = 


(0,0,P), 


p = 


(0,0,0), 


p = 


(p,p,o) 


p = 


(p,0,p) 


p = 


(0,p,p) , 



(25) 

with the indices of the operators chosen accordingly. 

The matrix elements of the operators given in Eqs. (|20|). (|22fl are then extracted from 

ratios of two-point functions. In forming the ratios from Eq. (|18|). we first need to determine tq 
for the various operators. We find Tn a = +, r n b = — , Tr> a = + and r n b = + 1731. 

41 ^44 412 ^411 

This gives the ratios (for <C t <C L t ) 

C°*t(t) 

R la = ^01^ = ^(0^ , (26) 

R2a = = ^ 2) " are ' (28) 

where i and j are spatial indices, and O4 = A^{x) = q(x)^f^^u{x) is the operator given in Eq. ()19|) 
with no derivatives and fx = 4. In Eq. (j2U), F(Ep, t) = coth [Ep(t - L t /2)\ and tanh [Ep{t - L t /2)] 
for J(x) = II(x) and J(x) = A±{x) pseudoscalar mesons, respectively. 

Figure n shows a typical example of the ratio in Eq. (ffify using a J(x) = A±(x) pseudoscalar 
meson ((£)a 5 )' where we observe a plateau between t = 7 and t = 40. After extracting R la from 
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FIG. 1: The imaginary part of R la as denned in Eq. l(2"fi|l using a J(x) = Ai{x) = q{x)^i"f^u{x) meson 
interpolating field, for (3 = 5.29, K sca = 0.13550 and valence masses, K va ii = 0.13550, K va i2 = 0.13430. 



the plateaus, we use Eq. (|26j) to extract (£)a are . Similarly, a hyperbolic tangent fit to the ratio 
R lb in Fig. [21 and a constant fit to R 2a in Fig. 01 allow for the extraction of (C)fe are an d (£ 2 )a are , 
respectively. 

Here and in the following, we use the notation (£ n ) 5 and (£ n ) 45 to distinguish the results for 
J(x) = n(x) and J(x) = A±(x) pseudoscalar mesons, respectively. 



D. Operator Renormalisation and Mixing 

In general, bare lattice operators must be renormalised in some scheme S and at a scale M. If 
the operator is multiplicatively renormalisable, which is the case for the operators (|20JI and (|21[) . 

we have 

O s (M 2 ) = Z$(M 2 )0(a) , (29) 

where 0(a) denotes the bare operator at lattice spacing a. Since (£ n ) is computed from a ratio of 
two-point functions with the operator O in the numerator and the 4-component of the axial vector 
current O4 in the denominator the renormalised value is given by 

{£ n ) S {M 2 ) = Z o( M2 ) ^bare ^ ^ 
Zo 4 



if O is multiplicatively renormalisable. 
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FIG. 2: R lb as defined in Eq. f2"7| using a J(x) = Ai{x) = q(x)j4j5u(x) meson interpolating field, for 
(3 = 5.29, K sca = 0.13550 and valence masses, K va ii = 0.13550, K va i2 = 0.13430. Fit function is y = 
Atanh[fe(i — L t /2)}, where A and b are fit parameters. 



TABLE II: Results for the non-perturbative RGI renormalisation constants, Z RG1 , for the operators defined 
in Eqs. 0, lEU and as well as for C 4 . 

a 7R.GI 7RGI 7RGI 7 

P Z la Z lb Z 2a Z 4 



5.20 1.52(4) 1.55(5) 2.4(1) 0.765(5) 

5.25 1.52(4) 1.55(5) 2.4(1) 0.769(4) 

5.29 1.54(4) 1.56(5) 2.45(10) 0.772(4) 

5.40 1.57(3) 1.60(4) 2.5(1) 0.783(4) 



In this work, we renormalise our operators non-perturbatively. Here we restrict ourselves to a 
short outline of the procedure. More details can be found in Section 5.2.3 of Ref. [64|. and a fuller 
account will be given in a forthcoming publication. 

We start from a MOM-like renormalisation condition imposed on the lattice jfjjj |(3f| and per- 
form a chiral extrapolation of the non-perturbative renormalisation factors at fixed [3 and fixed 
momentum. We then apply continuum perturbation theory to calculate the renormalisation group 
invariant renormalisation factor Z RGI from the chirally extrapolated Zs 64]. Our results for the 
operators and (|2T jl . i.e. Zf^ 1 and Zf b GI , can be found in Table [nj where also Zq 4 is given. 
Note that Z and Z RGl coincide for O4 because the anomalous dimension of the axial vector current 
vanishes. 



In the final step we have to convert Z RGI to the MS scheme at some renormalisation scale 
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FIG. 3: R 2a as defined in Eq. l(2*%|l using a J(x) = q{x)^§u(x) meson interpolating field, for j3 — 5.29, k s 
0.13590 and degenerate valence masses, K va ii = 0.13490, K va i2 = 0.13490. 



ji 2 . This is done perturbatively, and the result depends on the value of Atjc in physical 



M 2 



units. We use tqA- 



MS 



0.617 Vm and r = 0.467 fm to obtain A 



= 261 MeV. For the operators 
(HHP and dHJ) we find Z m /Z RGI = 0.7154 at the scale p 2 = 4GeV 2 . 

If there are operators having the same quantum numbers and the same or lower dimension, they 
may mix with the operator we are interested in and we must renormalise our operator via 



Of(M 2 ) = yzf 1 (M 2 ,a)0 J (a) 



(31) 



From |60J we know that 0% X2 
operator can be written as 



(Eq. (J22J)) mixes with O, 



9 (Eq. (|24|) ) such that the renormalised 



Q s 



412 



^412^412 



(32) 



If we restrict ourselves to forward matrix elements, the operator O ^^ cannot contribute and O412 
becomes effectively multiplicatively renormalisable. Thus we can compute Zf ?9 in Eq. (j3U) non- 
perturbatively as sketched above. A sample result is shown in Fig. 2 of Ref. |6J3] 



where Z^* 1 is 



called Z^qj- Our numbers for are also given in Table ITTI 

The mixing factor Z^ ix , on the other hand, has only been computed in one-loop tadpole- 
improved lattice perturbation theory jH^]. In order to avoid the logarithms in the perturbative 
expressions we work at the scale [i 2 = 1 /a 2 , where a is obtained from the value of ro/a in the chiral 
limit 63] . In this way we find the numbers given in Table II I II 
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TABLE III: Results for the renormalisation mixing coefficient, Z^^, computed in tadpole- improved pertur- 
bation theory in the MS scheme at the scale /x 2 = (1/a) 2 , where a is obtained from the value of ro/a in the 
chiral limit |6jj. 






H 2 = 1/a 2 [GeV 2 ] 


mix 


5.20 


5.3361 


-0.00258 


5.25 


6.2001 


-0.00253 


5.29 


6.9696 


-0.00250 


5.40 


9.7344 


-0.00240 



TABLE IV: Results for Z^ s /Zf^ 1 at /i 2 = (1/a) 2 for the lattice spacings a in our simulations. 

P — (1/a) 2 [GeV 2 ] zW/Zf? 1 
5.20 5.3361 0.5650 
5.25 6.2001 0.5545 
5.29 6.9696 0.5465 
5.40 9.7344 0.5262 



The values of the conversion factor Z^/Z^ 1 at the relevant scales are again computed in 
continuum perturbation theory and are collected in Table ITVl 



Denoting the unrenormalised values of /n and (£ 2 ) by f^ re and (£; 2 ) bare , respectively, we have 
from Eq. (|16j) 



(0\O2 l2 \U(p)) = /n">ip 2 p4<£ 2 > bare , (33) 



and 



(o\of 12 \n( P )) = f^ PmP4 (zf 12 <e 2 ) w + zt 

= fn P iP2 Pi (^(e) h -+^) . (34) 

Here the renormalised /n is given by 

fn = Z J^ e , (35) 

and for the renormalised (£ 2 ) we get 

(£ 2 ) = M2 (e 2 ) barc + ^ix > (36) 
ZQa ^04 

So we first obtain (£ 2 ) at the scale = (1/a) 2 . Using the relation between (£ 2 ) and the 
Gegenbauer moment a 2 , Eq. (fTT|) . along with the NLO scale dependence of the latter, Eq. ((7J), we 
get (£ 2 ) at the scale fj, 2 = 4 GeV 2 . We calculate the running coupling from the 4-loop approximation 
of the /3-function in the MS scheme with Aye = 0.261 GeV 



;oupii 
61]. 
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TABLE V: Bare results for (£ 2 };J and (£ 2 )a 5 calculated on each dataset with degenerate valence quark masses 



p 


^sca 




/ j-9\ 4^ 

(Z)a 


5.20 





13420 





1353(47) 


0.1447(46) 


5.20 





13500 





1296(77) 


0.1575(62) 


5.20 





13550 





1518(65) 


0.140(10) 


5.25 





13460 





1380(55) 


0.1328(82) 


5.25 





13520 





1450(67) 


0.1706(57) 


5.25 





13575 





1371(82) 


0.1541(93) 


5.29 





13400 





1434(54) 


0.1537(47) 


5.29 





13500 





1346(37) 


0.1587(35) 


5.29 





13550 





1578(76) 


0.1737(68) 


5.29 





13590 





1401(94) 


0.1769(77) 


5.40 





13500 





1488(42) 


0.1516(58) 


5.40 





13560 





1581(86) 


0.1780(74) 


5.40 





13610 





1495(83) 


0.172(11) 



III. NUMERICAL RESULTS 
A. Mesons with mass degenerate quarks 

Investigating quark mass degenerate mesons, i.e., the matrix element Eq. (|16|) using the operator 
in Eq. (|15j) with identical masses for the fermion propagators, allows us to investigate the structure 
of the pions. In this case, all odd moments vanish, hence we focus on the lowest non-trivial moment, 

For each of our datasets, we extract a value for (£ 2 ) bare from Eq. (|28|) and renormalise using 
Eq. Q36|). In TablelVlwe present our results for (£ 2 ) bare . "We find that the results for (£ 2 ) using the 
meson interpolating operator lead to very poorly constrained chiral and continuum extrapolations 
for operators involving 2 derivatives. Hence in the following we only discuss the results for (£ 2 ) 
obtained using the IT interpolating field. 

In order to obtain a result in the continuum limit at the physical pion mass, we first extrapolate 
our results at constant (3 to the physical pion mass. In Fig. 0] we display the chiral extrapolations 
for (3 = 5.40 (top) and 5.29 (bottom), while Fig. [21 contains the corresponding extrapolations for 
f3 = 5.25 (top) and 5.20 (bottom). These results exhibit only a mild dependence on the quark 
mass and their values at the physical pion mass agree within errors. The smooth linear behaviour 
of (£ 2 ) was predicted in Ref. pa . I&gj ] where it was shown that at next-to- leading order in chiral 
perturbation theory, all possible non-analytic corrections to the matrix elements (|lfi[) are contained 
in /n- 

Now that we have calculated results at the physical pion mass for each choice of /3, we are in 
a position to examine the behaviour of our results as a function of the lattice spacing. In Fig. El 
we use the values of ro/a extrapolated to the chiral limit for each (3 (see Table 3 of Ref. [6Jj|) to 
study the dependence of our results on the lattice spacing. Even though our operators are not 
C(a)-improved, we find a negligible dependence on the lattice spacing, at least when compared to 
the statistical errors. 

Employing a linear extrapolation to the continuum limit at the physical pion mass, we find the 
second moment of the pion's distribution amplitude to be 

(£ 2 >?V = 4 GeV 2 ) = 0.269(39) , (37) 
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FIG. 4: Chiral extrapolation of {£, 2 )-n at constant (3 for /3 = 5.40 (top) and /? = 5.29 (bottom) for Of 12 from 
Eq. J23 in the MS scheme at fi 2 = 4 GeV 2 . 



with an acceptable x 2 /dof = 0.5, which is close to the value (O^V = 4 GeV 2 ) = 0.286(49)1^ 
found in Ref. .,55], and larger than the asymptotic value, (£ 2 ) as = 0.2. 



B. Mesons with mass non-degenerate quarks 

When the masses of the quark and the antiquark in Eq. (|16[) become unequal, the odd moments 
will no longer vanish and — with appropriate adjustment of the quark masses — we can directly 
obtain the corresponding moments of the Kaon. The results that will be discussed in this section 
are tabulated in Tables I VI II IVIII1 IIXI and [Xj Due to the large amount of resources required to 
calculate these partially quenched results, we simulate at a fixed value of (3 = 5.29 where we 
have four different sea quark masses at our disposal. As a result, we are not able to examine the 
lattice spacing dependence of these results. However, we take encouragement from our results in 
the previous section, where we found that discretisation effects are small for (£ 2 ), and neglect the 
extrapolation to the continuum limit. 
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FIG. 5: Chiral extrapolation of (£ 2 )tt at constant (3 for — 5.25 (top) and (3 — 5.20 (bottom) for 0\ 12 from 
Eq. in the MS scheme at fi 2 = 4 GeV 2 . 



Occasionally the raw data is so noisy that it is not possible to perform a stable fit to one or 
more of the ratios in Eqs. (|26j). (f2"T|) and (|2S|) . In such instances, we are unable to report a result 
and hence gaps are present in Tables IV 111 IV1111 11X1 and [X] 



1. Second moment 

Figure [7| shows the second moment, (£, 2 )k, extracted from Eq. (|28|) at the working point, 
P = 5.29, ^sea — 0.13500, as a function of the squared Kaon mass, n-^, for various choices of the 
valence quark masses. Here when we refer to the Kaon mass, we mean a pseudoscalar mass which 
is a function of two valence quarks, tuk = w ps (K va n, ft va i2), where K va n > K va i2- (These masses 
are provided in the third columns of Tables IVIIfJXl ) A solid symbol indicates the point where 
K vaii = K vai2 = K sea- The behaviour towards the chiral limit of the available data points suggests 
that indeed a linear extrapolation is possible. The vertical dotted line indicates the physical kaon 
mass, rriK = 0.494 GeV. 
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FIG. 6: Results for (£ 2 )-n- for each value of (3 at the physical pion mass as a function of a 2 for 0% 12 from 
Eq. J23 in the MS scheme at fi 2 = 4GeV 2 . 
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FIG. 7: (£ 2 )if , extracted from Eq. (|2*%)) at the working point, (3 = 5.29, K sca = 0.13500, as a function of 
the squared Kaon mass, m 2 K , for various choices of the valence quark masses. Results are quoted in the MS 
scheme at /i 2 = 4 GeV 2 . The vertical dotted line corresponds to the physical Kaon mass. 



In order to obtain a result at the physical ir and K masses, we performed similar fits at all 
available sea quark masses corresponding to K sea = 0.13400,0.13500,0.13550,0.13590, and then 
tried to extrapolate in the sea quark mass (or = m ps (K seSu , n sea )) to = 0.140 GeV. This 
final extrapolation, however, turns out to be unreliable (large x 2 /dof). Hence we attempt to fit to 
all the data available with the global ansatz 



(£, 2 )k = a + /3m^(K sea , K sea ) + 7m| r («; V aii, K va i2) , 



(38) 
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FIG. 8: Results of fit in Eq. (|38|l for m n — m£ ysl , together with our results for (^ 2 )k from all four values 
of K sea at /3 = 5.29 considered. Results are quoted in the MS scheme at fi 2 = 4 GeV 2 . The vertical dotted 
line corresponds to the physical Kaon mass. 



with three fit parameters, a, /3, 7. After performing such a fit we find 

a = 0.264(7), (3 = -0.00005(841), 7 = -0.016(9) , (39) 

with a x 2 /dof = 1.06. The fit results indicate that the dependence of (£, 2 }k on the sea quark mass 
is negligible, while the dependence on the valence quarks is very small. 

In Fig. El we display all our results for (£, 2 }k for all four sea quark masses, together with the 
fitted ansatz, Eq. I|38j). at the physical pion mass, i.e. a + f3m 2 p j l y S +7Wi^-(/c va ii, fiVal2)> given by the 
solid line. For further clarification, the result of this fit is also shown in Fig. EJJfor each value of n sea 
separately. In this figure, each solid line corresponds to the fitted ansatz, Eq. (|38)l. for m 7r (K sea , K sea ) 
evaluated at, going from top to bottom, K sea = 0.13400, K sea = 0.13500, K sea = 0.13550 and 
^sea = 0.13590. For example, in the top figure, the solid line refers to a + (3m 2 n (0.13400, 0.13400) + 
7m| ( -(K va ii, «Val2); where m^(0. 13400, 0.13400) is taken from Tabled 

To obtain our final result, we insert the physical values for and mjf, together with the fitted 
parameters in Eq. (|3*9l . into Eq. PH)) and we find in the MS scheme at /x 2 = 4 GeV 2 

(£ 2 )IV = 4 GeV 2 ) = 0.260(6) . (40) 

Since we only have results with non-degenerate quark masses at one value of j3 = 5.29, we are not 
able to perform a continuum extrapolation of {£, 2 )k- We are, however, able to gain an estimate 
of the systematic error due to discretisation effects by comparing the result for (^ 2 ) 7r at f3 = 5.29 
with that in the continuum limit (|37j) . Such a comparison suggests that there is a systematic error 
of roughly 6% due to discretisation effects. 

Comparing the results in Eqs. (|37|) and l)40|) . we see that second moments for the Kaon and 
pion coincide within errors, in agreement with findings in Refs. jy], 3^|. 
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FIG. 9: Results for (£ 2 }k as in Fig. [5] but with all four values of K sea separated for clarity. Each solid 
line corresponds to the fitted ansatz, Eq. I|38|) . for ?n T (K soa , K soa ) evaluated at, going from top to bottom, 
K soa = 0.13400, K sca = 0.13500, K sca = 0.13550 and K soa = 0.13590. The vertical dotted lines correspond to 
the physical Kaon mass. 



18 








06 







05 







04 








A 





03 








V 











02 







01 
















06 : 







05 : 







04 : 








LO 

A 





03 : 








V 











02 : 







01 : 






: 



0.1 



0.1 



0.2 



0.3 



(m^-m 2 ) [GeV 2 ] 



0.2 



0.3 



(m^-m 2 ) [GeV 2 ] 



0.4 



0.5 




0.4 



0.5 



FIG. 10: Result for and (£}| for (3 = 5.29, K sca = 0.13500, in the MS scheme at = 4 GeV 2 . The 
vertical dotted line corresponds to the physical m 2 K — m 2 mass difference. 



2. First moment 

Fi gures ITT71 and 1111 show the first moment, for the working points j3 — 5.29, Ksea — 0.13500 

and = 5.29 , K sea — 0.13590, respectively, as obtained from the two operators (j2UJ) and O^a, 
(|21jl. The results are plotted as a function of the mass splitting of the two quarks making up the 
meson, or more specifically m 2 K — m 2 . Here mx is the mass of a pseudoscalar meson constructed 
with one heavy and one light quark, while m n is the mass of a pseudoscalar meson constructed 
with two light quarks, i.e. m^O^an, K val2 ), mJn ^ , K va u) with /t vall > K val2 . The points lie on 
a straight line, once again as predicted in Ref. [68L l69j] . 

The vertical lines in Figs. ^] and ^2 show the location of the physical K-tt mass splitting and 
it is here that we extract our results for {£)k at each sea quark mass. These results are given in 
Table IVD together with slopes obtained from the simple fit 



d) K = B{m 2 K -ml) 



(41) 
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FIG. 11: Result for (£) 5 a and for (3 = 5.29, K sca = 0.13590, in the MS scheme at = 4 GeV 2 . The 
vertical dotted line corresponds to the physical m 2 K — m 2 mass difference. 



We observe that at each K se& , the four sets of results obtained with two different operators and two 
different Kaon interpolating fields all agree well. 

In order to extract a result at the physical pion mass, we examine the sea quark mass dependence 
of our results by plotting them as a function of the pion mass calculated with K va i = K sea (Table IJ) 
in Figs. ^1 and for the operators O a ^ v (|2U)) and O ([21 j). respectively. We extrapolate linearly 
in the mass of the light quark to the physical pion mass and quote the results in the last row of 
Table EH Avera ging over the four results, we find 

(OWiv 2 = 4 GeV 2 ) = 0.0272(5) . (42) 

Similar to the result for {C 2 )k hi Eq- (|40[h we expect that there is a systematic error of roughly 
6% due to discretisation effects. 
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TABLE VI: (Ok a t the physical mj^ — m 2 . mass splitting together with the slope, B, of the fit in Eq. (1411) for 
(3 = 5.29, in the MS scheme at /i 2 = 4 GeV 2 . The last row contains the values for (Ok chirally extrapolated 
in the sea quark mass to the physical point. 
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0.13400 
0.13500 
0.13550 
0.13590 


0.0215(5) 
0.0234(2) 
0.0240(5) 
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0.0261(3) 
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0.0252(11) 
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FIG. 12: Sea quark mass dependence of for (3 — 5.29, in the MS scheme at fi 2 — 4 GeV 2 . The vertical 
dotted lines correspond to the physical tt mass. 



IV. SUMMARY AND CONCLUSIONS 

We have presented results for the second moment of the pion's distribution amplitide and the 
first two moments of the Kaon's distribution amplitude, calculated on lattices generated by the 
QCDSF/UKQCD collaboration with two flavours of dynamical fermions. We use nonperturbatively 
determined renormalization coefficients (apart from the mixing with the operators containing total 
derivatives, which is calculated perturbatively) to convert our result to the MS scheme at 4 GeV 2 . 
Our results give model-independent insights into the distribution amplitude of pseudoscalar mesons 
with degenerate and non-degenerate quark masses. 

We find for the pion (O 1 )-* = 0.269(39), which is in agreement with other results appearing in the 
literature and larger than the asymptotic value. For the isT-meson we obtain (Ok = 0.0272(5)(17) 
and ({, 2 ) K = 0.260(6) (16), where the first error is statistical and the second is an estimate of the 
systematic error due to the fact that we have results with non-degenerate quarks at one value of 
P = 5.29 only, i.e., no continuum extrapolation. 

The coefficients a n in the Gegenbauer expansion of the DAs in Eq. (J3J) are related to the moments 
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FIG. 13: Sea quark mass dependence of for /3 = 5.29, in the MS scheme at /i 2 = 4 GeV 2 . The vertical 
dotted lines correspond to the physical tt mass. 



(£ n ) by simple agebraic relations (|17jl. Using our result in Eq. (j37j) we obtain, for the tt- meson 

a£(^ 2 = 4 GeV 2 ) = 0.201(114) , (43) 

and from Eqs. (|40j) and (|42|) for the i^-meson: 

of Gu 2 = 4 GeV 2 ) = 0.0453(9) (29) , (44) 
of (/u 2 = 4 GeV 2 ) = 0.175(18)(47) . (45) 

While our result for a\ is larger than the transverse lattice result |42|], all three numbers are 
well within the range suggested by QCD sum rule estimates and supported (for the pion) by the 
analysis of CLEO data on the ^7*7 transition form factor, cf. Eqs. (|llj) and (|12j) . Also the SU[3) 
breaking in the second Gegenbauer coefficient turns out to be small, in agreement with pul . yj). 
We note that in the context of SU(3) flavour violation, one might be worried about the absence of 
a dynamical strange quark in our simulations, however there has recently appeared a Nf = 2 + 1 
lattice calculation of 70] which is in good agreement with our result, giving us confidence that 
the effects of a dynamical strange quark are probably small. 

Our results indicate that it is important to consider not only the chiral extrapolation of the 
lattice results to the physical quark masses, but also to perform simulations at small enough lattice 
spacings to allow for a reliable extrapolation to the continuum limit. 

The corresponding DAs obtained by the truncation of the general expression in Eq. © after 
the second term are shown in Fig. Eland Fig. El for the tt and the iT-mesons, respectively. Note 
that the fT-meson DA is tilted towards larger momentum fractions carried by the heavier strange 
quark, which is in agreement with general expectations. 

In order to illustrate the possible effect of higher-order terms in the Gegenbauer expansion, we 
also show in Fig. El the pion DA obtained with the addition of the fourth-order polynomial with 
the coefficient a\ = —0.10(5) taken from Ref. jj^]. In both cases (with and without a\) the value 
of the DA in the middle point agrees well with the estimate in Eq. (|13|) . The question whether the 
"camel-hump" structure of the DA is present in the physical DA depends on the contribution of 
yet higher-order polynomials that are beyond the reach of the present analysis. 
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FIG. 14: Distribution amplitude of the pion using the expansion in Eq. with our result for a\ — 0.201(114) 
and a\ = 0. This result is obtained in the MS scheme at \j? — 4GeV 2 . The shaded area indicates the results 
obtained when varies between the maximum and minimum values allowed by its error. 
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APPENDIX A: LATTICE RESULTS BY WORKING POINT 

The following tables summarize our findings at individual working points. 
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FIG. 16: Distribution amplitude of the pion using the expansion in Eq. J2J with = 0.201(114) and 
a\ = —0.10(5). This result is obtained in the MS scheme at fi 2 = 4GeV 2 . The shaded area indicates the 
results obtained when we fix a\ at its central value and let aj vary between the maximum and minimum 
values allowed by its error. 
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TABLE VII: Bare results for (£) and (£ 2 ) for various choices of valences quarks with K va n > K V ai2, together 
with the corresponding pseudoscalar mass, m ps (K va n, K V ai2)- Here (3 = 5.29 and K sca = 0.13400. 
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0.13400 



tot 



^val2 
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0.13580 
0.13630 
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0.4541(11 
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0.3913(14 
0.4639(15 
0.4423(16 
0.3929(21 
0.3569(32 
0.305(14) 



-0.00027(99) 
0.0053(10 
0.0179(14 
0.0121(14 

0.0259(15 
0.0201(16 
0.0076(19 

0.0332(19 
0.0275(19 
0.0154(23 
0.0081(28 



0.0001(11) 
0.0059(11) 
0.0190(15) 
0.0131(15) 
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0.1607(70) 
0.1621(75) 
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0.1633(91) 
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TABLE VIII: Bare results for (£) and (£ 2 ) for various choices of valences quarks with K va n > K va i2, together 
with the corresponding pseudoscalar mass, m ps (K va n, K va i2). Here [3 = 5.29 and K sca = 0.13500. 
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TABLE IX: Bare results for (£) and (£ 2 ) for various choices of valences quarks with K V aii > ^vai2, together 
with the corresponding pseudoscalar mass, m ps (K V aii, « va i2). Here (3 = 5.29 and K sca = 0.13550. 



= 



5.29, K sca = 0.13550 

(Of 



<or 



(0T 



0.13430 
0.13430 
0.13490 
0.13430 
0.13490 
0.13550 
0.13430 
0.13490 
0.13550 
0.13585 
0.13430 
0.13490 
0.13550 
0.13585 
0.13620 



0.46159(55) 
0.43065(56) 
0.39820(58) 
0.39837(59) 
0.36397(60) 
0.32696(64) 
0.37873(62) 
0.34288(64) 
0.30381(67) 
0.27895(72) 
0.35835(74) 
0.32065(78) 
0.27880(88) 
0.2515(11) 



-0.0006(8) 
0.0109(9) 
-0.0008(11) 
0.0226(11) 
0.0106(14) 
-0.0020(16) 
0.0293(13) 
0.0166(18) 
0.0045(19) 
-0.0034(22) 
0.0355(17) 
0.0220(24) 
0.0098(26) 
0.0038(27) 



0.2198(26) -0.0081(56) 



-0.0005(10) 
0.0110(12) 
-0.0008(14) 
0.0228(14) 
0.0106(17) 
-0.0019(20) 
0.0296(16) 
0.0166(21) 
0.0049(23) 
-0.0026(27) 
0.0362(21) 
0.0226(27) 
0.0112(30) 
0.0051(33) 
-0.0044(51) 



0.0115(1) 0.0116(2) 



0.0234(4) 
0.0118(2) 

0.0304(6) 
0.0190(5) 
0.0071(3) 

0.0380(13) 
0.0266(13) 
0.0147(12) 



0.0236(4) 
0.0120(3) 

0.0306(8) 
0.0190(7) 
0.0071(6) 

0.0382(17) 
0.0269(21) 



0.1484(37) 
0.1510(41) 
0.1525(45) 
0.1546(47) 
0.1559(52) 
0.1578(76) 
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0.162(10) 
0.1613(64) 
0.1646(88) 
0.167(11) 
0.168(13) 



0.1579(37) 
0.1617(40) 
0.1645(45) 
0.1670(46) 
0.1698(52) 
0.1737(68) 
0.1713(53) 
0.1739(66) 
0.1791(78) 
0.1841(91) 
0.1770(66) 
0.1812(83) 
0.1870(97) 
0.193(11) 



0.173(33) 0.201(18) 



TABLE X: Bare results for (£) and (£ 2 ) for various choices of valences quarks with K va n > K va i2, together 
with the corresponding pseudoscalar mass, m ps (K va ii, K va i2). Here [3 = 5.29 and K sca = 0.13590. 



5.29, Ksca 

(0^~ 



0.13590 



(or 



0.13490 
0.13490 
0.13530 
0.13490 
0.13530 
0.13575 
0.13490 
0.13530 
0.13575 
0.13590 
0.13490 
0.13530 
0.13575 
0.13590 
0.13617 



0.37239(65 
0.34894(67 
0.32430(69 
0.32121(71 
0.29483(73 
0.26270(77 
0.31158(74 
0.28450(75 
0.25123(80 
0.23925(82 
0.29364(80 
0.26503(81 
0.22921(86 
0.21603(89 
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-0.0018(20) 
0.0197(13) 
0.0110(15) 
0.0008(19) 
-0.0024(25) 
0.0244(17) 
0.0157(19) 
0.0054(25) 
0.0019(28) 
-0.0042(44) 
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0.0062(11) 
-0.0029(14) 
0.0156(13) 
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0.0186(15) 
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-0.0042(26) 
0.0238(20) 
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